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SIMULATION OF A TURING MACHINE 
ON A DIGITAL COMPUTER 


Robert W. Coffin,* Harry E. Goheen,** and Walter R. Stahl*** 


I. INTRODUCTION 


The theory of algorithms relies heavily upon 
the conceptual and theoretical usefulness of the 
Turing Machine.“ Recent work by Trakhen- 
brot (? has given further support to the hy- 
pothesis that “all algorithms can be given in 
the form of functional matrices and executed 
by the corresponding Turing Machines.” Such 
a statement does not immediately suggest that 
all problems should be reduced to their equiva- 
lent Turing Machine, but the implication is 
clear that if certain problems, recognizable as 
algorithms, do not lend themselves to a solution 
in a formal logic structure, they may be reduced 
to a Turing scheme using a suitable strategy. 
Examples of such problems are revealed in the 
work by Lusted and Stahl) in the problem of 
medical diagnosis and by Stahl and Goheen“? 
in simulation of the operation of biological cell 
systems. 


For years the practical value of the Turing 
Machine has been discounted because there is 
widespread belief among mathematicians that 
a Turing Machine always requires a tape of 
infinite length and that there is no assurance 
that a particular algorithm will achieve a stable 
solution in a finite length of time. These mis- 
conceptions probably stem from the fact that 
the Turing Machine invariably enters the dis- 
cussion of complex, self-organizing automata 


in which these conditions might exist. Further- 
more, the work of A. M. Turing’ preceded the 
growth of modern electronic technology by 
several years, consequently to implement a 
Turing Machine as a physical device, at that 
time, would have been an impractical and costly 
undertaking. 


The authors and their colleagues have found 
that simulation of the Turing Machine on a 
Digital Computer is a useful and practical tool 
not only for problem solving and validation of 
algorithms but also for teaching students the 
fundamentals of programming. This paper will 
describe the digital program by which a gen- 
eralized simulation of the Turing Machine has 
been accomplished. 


II. PRELIMINARY CONSIDERATIONS 


If there are M configurations corresponding 
to the Q-levels of the Turing Machine each re- 
quiring m bits, and if the Turing Machine 
requires an alphabet of size N, each symbol 
requiring n bits, we shall need a list of MN 
strings which we call machine quintuples. A 
machine quintuple will consist of : 

(a) m bits identifying a Q-level, 

(b) n bits identifying a symbol of the 
alphabet, 

(c) m bits designating the next Q-level, 
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(d) 2 bits designating the next motion 
(either Left, Right, or Place), 


(e) n bits designating the symbol to 
replace the symbol represented by (b) 
at the current location in the Turing 
tape. 


In the descriptive language of the Turing 
Table, the last (m + 2 + n) bits represent the 
(i, j) entry of the Turing Table, where i and j 
are given in (a) and (b) of the quintuple. The 
list of quintuples will require a maximum of 
2MN (m+ n + 1) bits. As a practical matter, 
Turing Tables are rarely “saturated,” i.e., the 
entire matrix filled with non-zero entries, how- 
ever the maximum bit requirement: 


2MN (m+n+4+1)=A (1) 
is an important program design criterion. 


To provide for rapid and intelligible inter- 
pretation of a processed Turing tape, it is 
desirable to limit the magnitude of the Turing 
alphabet to those alphabets of available input- 
output devices. It is important to note that this 
restriction is one of convenience since, in binary 
form, a symbol is limited only by the allow- 
able n. If a Turing Machine requires a large 
alphabet it would be necessary to numerically 
precode the alphabet, or alternatively, use a 
two-for-one notation. 


In addition to (1) a portion of the total 
available memory must be used to store the 
Turing tape. If this memory allocation is b 
bits, then the maximum possible length of the 
Turing tape will be the parameter: 


[b/n] = B (2) 


If the total number of bits in the memory of 
a specific computer is S, and the simulation 
program requires s bits, we have the following 
inequality which defines the maximum com- 
bination of the Turing Machine and tape that 
can be processed. 
A+B<S—s (3) 
In addition to (3), there is another restric- 
tion imposed upon the simulation program. We 
denote by r the time used in passing from one 
quintuple to the next. If an algorithm requires 
T quintuple cycles for solution, we have the fol- 
lowing inequality in which C denotes the time 
available for the processing of the algorithm. 


far 


TSC (4) 


The inequalities (3) and (4) are viewed as 
general restrictions upon the magnitude of 
Turing Machine which can be simulated on a 
particular computer. 


III. STRATEGY OF SIMULATION 


The simulation of a machine in a digital com- 
puter requires two major routines that are 
logically distinct from each other. The first 
of these routines serves to receive and interpret 
external notation and generate the “compiled 
logic” of the device to be simulated. Whether or 
not this routine is a true compiler, in the sense 
of current usage of the word, depends upon the 
logical complexity of the simulated device and 
the corresponding notation required to fully 
describe its function. The logic of the Turing 
Machine is quite simple and completely em- 
bodied in the construction of the quintuple 
which in turn is determined by the author of 
the Turing Machine. Therefore the builder 
is required only to make the notational transla- 
tion of externally coded quintuple statements 
into machine language. To distinguish between 
a true compiler and the routine that generates 
the “compiled logic”, we shall henceforth refer 
to the latter as the “builder” routine. 


The essential steps of the builder routine are 
outlined below. 
1. Read quintuple card. 


2. HALT card? Yes, goto 14. 
3. STOP quintuple? Yes, go to 12. 
4. Extract (Q, S, Q, M, S). 
5. Convert Q-s to binary. 
6. Assign movement code. 
7. Form function: (Q, M, S). 
8. Form identifier (Q, S). 
9. Store identifier and function. 
10. Update storage. 
11. Go to 1. 
12. Insert stop code for function. 
13. Go to 8. 


14. Exit Builder. 


Each quintuplet is punched on a separate card 
in the format 
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where 10 A is the Turing Matrix identification 
pair and 11 R B is the function triplet. 


The only exception to this format is the 
special stop card with the format: 


QSSTOP 


We have chosen to use a special stop code 
for three reasons. First, the special stop code 
makes it possible to distinglish between a nor- 
mal halt and an error in the Turing Machine. 
Secondly, it is far more economical of computer 
time to make a simple, one-step test rather than 
go through three, multi-step tests to determine 
the (Qi, Sj, Qi, P, Sj) quintuple. Finally, the 
special stop code and the natural stop quintuple 
are theoretically equivalent. 


The quintuple deck is terminated with the 
halt card: 
HALT 


The dimensions of the Turing Table are de- 
termined by the builder. The row dimension is 
equal to the largest defined value of Q, and the 
column dimension is equal to the number of 


different symbols defined in the quintuple list. 


It must be noted that the quintuple notation 
makes the determination of the dimensions un- 
necessary except for purposes of arraying the 
matrix on a line printer or other device. 


Since the speed of the builder routine is lim- 
ited by the speed of the card reader there is 
time for nonessential checking for mispunched 
cards. The inclusion of three such checking 
features has proved to be very helpful in ridding 
quintuple decks of gross errors. The first and 
most valuable check is dependent upon an al- 
phabet definition card that is read immediately 
prior to the first quintuple card. A short routine 
generates a check list consisting of all of the 
defined symbols which can legally occur in the 
quintuple deck. Then as the quintuple cards are 
read, each symbol is checked against the defined 
alphabet. In the event an undefined symbol 
appears, the illegal symbol and the card contain- 
ing the error are typed out with an error mes- 
sage. The error routine keeps a list of the ad- 
dresses where the incorrect cards should be 
located, and the builder proceeds to the next 
card. After the HALT card is detected the pro- 
gram stops and allows the operator to insert 


the corrected cards in the card hopper. Upon 
restarting, the program will translate the cor- 
rected cards and store the corresponding entries 
in their original order. Similar checks are made 
to detect illegal movement and non-numeric 
symbols in the Q state, with identical proce- 
dures for error recovery. 


After the building is complete, before any 
simulation can take place, a Turing tape must 
be loaded into memory. The Turing tape is 
punched on one or more cards, and each symbol 
is checked against the original alphabet list. 
If an illegal symbol occurs in the Turing tape, 
an error message is typed and the reader 
routine halts. Upon restarting, the corrected 
tape is read, checked and stored. 


The second part of the program is called the 
“driver”? routine which operates on an element 
of compiled logic, provides functional continuity 
between elements, and hence forces the simula- 
tion. In the ideal case, no part of the driving 
routine performs any function except to operate 
on an element or provide continuity between 
elements, however, the timing restriction will 
certainly have to be included in the driver 
routine and is considered to be an allowable 
artifact. Other artifacts may be justified for 
purposes of demonstration or instruction but 
are undesirable burdens on the simulation since 
time is at a premium and each artifact must 
consume time on every quintuple cycle. Clearly, 
any artifact that causes the driver routine to 
amend itself or alter an element of compiled 
logic is unallowable. 


The following outline shows the essential 
steps of the driver routine. 


1. Extract symbol from Turing tape. 

2. Merge symbol with previous Q-level (from 
step 7). 

3. Find function corresponding to (Q, S). 

4. IsfunctionaSTOP? Yes, go to 14. 

5. Substitute function symbol in Turing tape. 

6. Move Turing tape “window”. 

7. Extract next Q-level. 

8. Has time been exceeded? Yes, go to 10. 

9 


Go to 1. 
10. Type time warning. 
11. Halt. 
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12. Insert new time restriction. 


13. Gotol. 
14. Output final Turing tape. 
15. Halt. 


The step which consumes the largest segment 
of time in the “Driver” is the table lookup to 
find the next function triplet, (step 3). The 
time required for this search can be greatly 
reduced if we take advantage of the natural 
tendency of Turing Machine authors to group 
important Q-levels together. Thus, the prob- 
ability of finding the next quintuple in a close 
proximity to the current quintuple is much 
higher than finding it at a great distance. After 
a check is made against the current quintuple, a 
flutter search is used in which the search begins 
with the next entry forward of the current 
quintuple and then shifts to the last entry 
preceding then to the second entry forward and 
so forth until the required quintuple is located. 
It is important to note that the order in which 
the quintuples occur is important only to the 
efficiency of the algorithm but does not influence 
its logical operation. 


IV. SAMPLE ALGORITHMS 


The first problem is to detect the change of 
pattern in a strip of contrasting white and black 
segments as shown in Fig. 1(a). 


Figure 1 (a). 


— . . 
AAAABBBAAABBAAABBB 


Figure 1(b). 


L WW EN E) 
AAAABAABAABAB AABAA 
Figure l(c). 


To code the pattern as a string of symbols we 
simply assign the letter A to the white segments 
and the letter B to the black segments, and in 
both cases make the number of letters propor- 
tional to the length of each segment as shown in 
Fig. 1(b). The problem is to write an algorithm 
in Turing notation that will operate on this 
coded string in such a way that only the points 
at which the pattern changes are marked. The 
desired solution is shown in Fig. 1 (c). This is 
logically equivalent to generating a unit pulse 
based upon a criterion of wave height or 
frequency. 


The Turing Table for this algorithm is shown 
in Fig. 2, and is a simple 3 X 3 matrix with 2 
null entries. 


k A B 
a Ra 
Q2 |2Px*/2RA|3RB 
as 


Figure 2. 


The strategy is to proceed from the left * to the 
right passing from Q1 to Q2 when either an A 
or B is encountered. If an A is seen, control 
stops in Q2 until a B is seen. The first B is 
retained and control goes to Q3 where succes- 
sive B’s are replaced by A’s until an A is seen. 
The first A is replaced by a B and control passes 
back to Q2 which begins the algorithm again. 
If a B is seen as the first symbol, the same type 
of procedure is followed except control goes im- 
mediately to state Q3. 


Fig. 3 shows the complete input deck as it 
appears before the start of a run. 
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*AB 

| * 2RA 

2 * STOP 

2A 2 R A 

2 8 3 R B 

3 * STOP 

3 A 2R B 

3B 3 R A 

HALT 

*AAAABBBAAABBAAABBB X 

END TEST 
Figure 3. 


The first card is the alphabet definition card 
consisting of the three legal symbols: *, A, and 
B. Following the alphabet card is the seven- 
card Turing program which is terminated with 
a HALT card. The coded test string shown in 
Fig. 1(b) is next, followed by an END OF 
TEST card. 


Fig. 4 shows the sequential operation of the 
Turing Machine. The current quintuple is 


*AAAABBBAAABBAAABBB* (661,*] 1892, R. *] 
192, Al 1992, R. Al] 
192, A] (992, R, Al 


1692, AJ 1862, R. A] 


*AAAABBBAAABBAAABBB* 
+AAAABBBAAABBAAABBB* 
+AAAABBBAAABBAAABBB* 
1692, A] 1662, R. Al 
1092,81 L683, R., B] 
1663,81 1663, R. A! 
1963,81 1863, R., A] 


+AAAABBBAAABBAAABBB* 
+AAAABBGAAABBAAABBE* 
+AAAABABAAABBAAABEB# 
*AAAABAAAAABBAAABBB* 
A [ØØ3,A] [ØØ2,R,B] 
(ØØ2,A] [ØØ2,R,A] 
1992, A] 1902, R, A] 
1692,81 1863, R, 8] 
1693,81 1603, R. A] 


*AAAABAABAABBAAABBB* 
*AAAABAABAABBAAABBB* 
+AAAABAABAABBAAABBB* 
#AAAABAABAABBAAABBB* 
+ AAAABAABAABAAAABBB* 
+ AAAABAABAABABAABBB* (ØØ3,A] 1682, R. 8] 
102, A] 1062, R, Al 
1882, AJ 1682, R. Al] 
1692,81 (083, R. 6] 
1603.81 (063, R. A] 
1693,81 1083, R. A] 
1663. *] STOP 


A 
*AAAABAABAABABAABBB* 
A 
*AAAABAABAABABAABBB* 


+AAAABAABAABABAABBB* 
+AAAABAABAABABAABAB* 
+ ARAABAABAABABAABAA* 
*AAAABAABAABABAABAA* 


RUN COMPLETE 
Figure 4. 


shown at the right side of the tape. The delta 
symbol under each line shows the next symbol 
to be viewed by the “window”. 


In this example the Turing tape is printed at 
the completion of every quintuple cycle con- 
sequently the “window” is shown to move only 
one symbol per line of output. Such an exten- 
sive output is desirable only for short, illustra- 
tive examples since an algorithm may involve 
many thousands of cycles for solution. 


The second algorithm is a popular child’s 
problem in which a number of boys and soldiers 
are to be transported across a river in a row 
boat. The boat is capable of holding one or two 
boys, or one soldier, but not one boy and one 
soldier. The basic solution strategy is a simple 
five step algorithm as follows: 

1. two boys across river 

Z. one boy back 

3. boy out, soldier in, soldier across. 
4. boy back 

5. gotol. 


When all of the soldiers are across the river, 
the algorithm reduces to a back and forth shut- 
tle involving only boys with a net gain of one 
boy per two river crossings. The problem can be 
coded in an interesting manner with symbols 
representing the various objects involved. 

G isa grassy spot on the river bank, 

W isthe water, 

F isa seat in the row boat, 

B isa boy, 

S isasoldier, 

% is the left bank of the river, 

} is the right bank of the river 

Æ is the tape end marker symbol 

— is an intermediate symbol used to 

denote movement of the boat. 


A sample configuration is shown in Fig. 5(a) 
with the solution configuration in Fig. 5(b). 


At the beginning all of the boys and soldiers 
are on the right side of the river with the boat 
at the right bank. In the final string, all of the 
boys and soldiers have been transported to the 
left side and the boat is at the left bank. 


Figure (6) shows the Turing program in 
quintuple form required to process any such 
input configuration. 
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#GGGGGX%WWFFOSBBSSG# 


Figure 5(a). ` 


* 388 S8 f FWWI 6666664 


Figure 5b). 
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Figure 6. 


2 IRZ 
3 3RF 
A # STOP 

4 SRG 
5 SRG 
6 6LS 
6 6L# 
7 7LB 
8 8RW 
8 8 10 L. 
9 8 LF 
9 8LB 
10 - 10 L 8 
10 B 
11 F 
12 * 
12 


ra = 
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In addition to the standard solution al- 
gorithm, all contingencies are anticipated such 
as the presence of only one boy; or the absence 
of soldiers. 


Because of the number of quintuple cycles 
involved it is not practical to show the test 
string after each cycle, however, figure 7 shows 
every second cycle. 


Note that the “window” appears to take a two- 
symbol jump between lines. 


V. DISCUSSION 


The ,sample algorithms were run on a ma- 
chine with a cycle time of 8 microseconds 
(SDS 920). The limiting factor, as mentioned 
above, is the time involved in finding the next 
quintuple in the total list. Since the efficiency is 
influenced by the order in which the quintuple 
deck is organized, the question arises whether 
or not it is possible to obtain an optimum order- 
ing. Generally, this would be quite unlikely 
since an optimum order is dependent on the 
input string which may be composed of a 
random group of symbols. However, if quin- 
tuples, containing symbols which are used in 
conjunction with each other, are placed to- 
gether in the input deck, the efficiency can prob- 
ably be improved. 


In the Boys and Soldiers algorithm we 
achieved a speed of 3,000 quintuple cycles/sec- 
ond by ordering the deck row-by-row as it 
appears in a Turing Functional Matrix. By 
rearranging some of the cards, particularly in 
Q-levels which deal with the left and right 
movement of the boat, we were able to increase 
the speed to about 3,300 quintuple cycles/sec- 
onds. In the opposite direction, we tried to 
arrange the deck in such a way that we pro- 
duced a very inefficient operation and were able 
to reduce the speed to slightly less than 1,000 
quintuple cycles/second. By shuffling the cards 
like a bridge deck, we obtained about 2,100 
quintuple cycles/second. 


It is not surprising that the row-by-row 
ordering is quite efficient since a Turing pro- 
grammer normally adopts a row-by-row strat- 
egy in solving a problem, consequently the 
probability is high that significant quintuples 
will be grouped together. Experience with a 


#GGGWZIWFFSBSBE 
UU FSS B= 
„add vf FBS B- 
SGCOUHUFFTRSBE 
ct Lesbo 
IGGGINWFFEGSB? 
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ICGOIHFFIGSGI 
IGGGZHWBBICSGI 
ecacmuitest:esce 
sqeoms-snesce 
PGCGINBBWICSCE 
rocowBBWcscE 
e 
IGEGxE-HWCSGF 
PE 
IGOGSFBNWEGSG? 
IGGBEFBHNNCSGY 
ab Eins 
„ach f BET Sa= 
WER ER Eser 
PE 
IOGRAFRUEGSCY 


A 
660 AF eb GC? 
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1801.6] 
1602.0] 
1063. 1 
1803, f] 
1003. L! 
1065.8] 
1006,81 
1606, 1 
1007. f] 
1008, 8! 
1912.—1 
(909,-) 
loss, u] 
1009, NI 
1912,38] 
1009, 8] 
1913,68] 
101,6] 
1016, 
1016,06 
1020, hl 
(022,F] 
1021.8] 
(019,58) 


Figure 7.1. 


(062,R,C) 
(602,R,C) 
(0935,R,W) 
(063,R,F) 
Loo, R. LI 
1005, R, S] 
1006.1. 8] 
1006, 1. L! 
1007. L. 61 
1099.1. 
1612. R. 1 
1909.1, 8] 
1908, R. u] 
1912, R. B] 
1099, , v1 
1908, 1, 8] 
1010, , f! 
1016, R. 3] 
1016, f. F! 
1019, 1, u] 
1022.1, 3] 
1021. 6, %! 
1019, R. E] 
1020, E,] 
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+GSGMINEF BGC 
vas FSG 
ash FoB80= 
76s F BG0· 
5s FEGBGO 
esch EG- 
ash HO d. 
roscuasriccas 
*cSGmiB-sriccce 
LES MES 
IGSGIWBEWIGCGY 
IGSGH-BWIIGGG? 
EREN C- 
LESSER: 
PGSGRBBWNTCGGE 
5s Cf nns 
„as Bf BHE. 
5asB F BOC 
„SB AF-HHu gad 


ray 
7685 f- BH .- 
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(619, f] 
(622,-) 
ETA 
(619, 4) 
(923,8) 
[024,#) 
182, 61 
1925, 61 
100,8] 
166,6] 
1606, f! 
1607. * 
1069, 
1012,81 
1899.8] 
loss, B] 
1912.—1 
(669,-] 
1968,71 
19010, 
(616, %) 
(016,8) 
(019,8) 
[922,-) 


Figure 7.2. 
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(022,-] 
(021,°) 
1919, ] 
192,111] 
1925,11] 
1026,81 
1927.7] 
1026, f] 
Loos, f] 
19012. 
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1020, f] 


Figure 7.3. 
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EGSR oh FIB C- 10 30, 11 
e 1930, f] 
eT rer 1830,81 
USRR CE [028,F] 
Ree 1029, . 
ee {611,W) 
VEER WE MOER (012,8) 
IGSBAWFBWIGGG? 1069, EI 
ENEE 1008, f] 
Keen 1012,-] 
ray 
Figure 7.4. 


(621,R,F) 
(619,L,W] 
1922, L, f! 
1821, R. W) 
1919, R. B] 
1826, R. F] 
1625, L.] 
1926, R. B] 
1926, R., 6] 
1928, 1, 61 
1929, L, B] 
1931. L, f] 
1931. L, u 
1030, R, 7] 
10 39, R. M1 
1630, R. f] 
1028, L,G) 
1629, L, B] 
1008, R. u! 
1012, R. F] 
1009, L. W] 
1608, 1, f] 
{611,L,-] 
1012, R.] 
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#GSBEFBWHIGGG# 1409, -1 1009, 1, 8] 
IGSBEFBUWEGGGY 1008,71 [Ø13,R,#] 
sasBerEWUZceGe 1013,81 [014,L,F) 
#GSBEFFWNIGCOF (014,#] (631,L,%) 
IGSBLFFUNZGCF 1031,81 (631,L,S] 
4885 Fe GGO- 1032,51 STOP 
RUN COMPLETE 
Figure 7.5. 


particular algorithm would certainly indicate 
an efficient input order. However, it is an in- 
teresting characteristic of a Turing program 
that it will run to solution regardless of the 
input order of the instructions. The authors 
know of no other programming language in 
which this is true. 


VI. CONCLUSION 


The value of the simulation of a Turing Ma- 
chine is two-fold. First, the simulation repre- 
sents a practical means of solving algorithms 
in their most fundamental form, that is, in a 
basic language which is independent of both 
computer and computer language. Secondly, 
for non-numeric problems the processing rates 
are found to be competitive with the more con- 
ventional computer languages. 
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